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ABSTRACT 

We investigate the role of the channel on H 2 molecule formation during the col¬ 
lapse of primordial gas clouds immersed in strong radiation fields which are assumed 
to have the shape of a diluted black-body spectra with temperature T^ad- Since the 
photodissociation rate of depends on its level population, we take full account of 
the vibrationally-resolved kinetics. We find that in clouds under soft but intense 
radiation fields with spectral temperature T^ad 7000 K, the channel is the domi¬ 
nant H 2 formation process. On the other hand, for harder spectra with Tj-ad ^ 7000 K, 
the H“ channel takes over in the production of molecular hydrogen. We calculate 
the critical radiation intensity needed for supermassive star formation by direct col¬ 
lapse and examine its dependence on the level population. Under the assumption 
of local thermodynamic equilibrium (LTE) level population, the critical intensity is 
underestimated by a factor of a few for soft spectra with Trad 7000 K. For harder 
spectra, the value of the critical intensity is not affected by the level population of H J. 
This result justifies previous estimates of the critical intensity assuming LTE popula¬ 
tions since radiation sources like young and/or metal-poor galaxies are predicted to 
have rather hard spectra. 

Key words: quasars: supermassive black holes-molecular processes-cosmology: the¬ 
ory. 


1 INTRODUCTION 

The collapse of primordial gas clouds can lead to the for¬ 
mation of ordinarily-massive stars (10 — 100 Mq, hereafter 
Pop III stars; see, e.g., Yoshida et al. 2008; Hosokawa et al. 
2011) or supermassive stars (M > 10® Mq, hereafter SMS; 
see, e.g., Bromm & Loeb 2003), depending on the strength 
of the radiation fields in the ambient medium. Pop III stars 
form by H 2 cooling in the absence of a strong Lyman-Werner 
(LW) radiation field {hu = 11.2 — 13.6eV). Conversely, in a 
cloud irradiated by an extremely strong LW radiation field, 
H 2 is completely destroyed, and cooling is suppressed. In 
the latter case, if the cloud is in a halo with virial temper¬ 
ature higher than ~ 10^ K, it can still collapse via atomic 
cooling (Omukai 2001), without experiencing major episodes 
of fragmentation (Bromm & Loeb 2003; Regan & Haehnelt 
2009a,b; Inayoshi et al. 2014; Latif et al. 2014; Regan et al. 
2014a,b; Choi et al. 2015; Becerra et al. 2015), and an em- 
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bryonic protostar can grow rapidly to a SMS by subsequent 
accretion (Hosokawa et al. 2012, 2013; Sakurai et al. 2015; 
Shlosman et al. 2015). The latter eventually collapses due 
to the post-Newtonian instability, leaving a black hole with 
~ 10 ®Mq (see, e.g., Shapiro & Teukolsky 1983). This pro¬ 
cess, called “direct”, or “monolithic” collapse, leads to the 
formation of SMS remnants that are promising candidates 
for being the seeds of the supermassive (> 10® Mq) black 
holes (SMBHs) observed at redshift 2 > 6 (see, e.g.. Fan 
et al. 2001; Mortlock et al. 2011; Venemans et al. 2013; Wu 
et al. 2015)). 

The feasibility of this scenario can be tested by com¬ 
paring the observed number density of the high -2 SMBHs 
n-SMBH with the theoretical predictions. By considering the 
probability distribution of the LW intensity J 21 around ha¬ 
los^ , nsMBH can be estimated from the critical value of LW 
intensity J 21 needed for the formation of SMS (Dijkstra 


^ We use the specific intensity at the center of the LW bands. 
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et al. 2008; Agarwal et al. 2012; Dijkstra et al. 2014; In- 
ayoshi & Tanaka 2015). Sugimura et al. (2014) (hereafter 
SOI14) found J 21 ~ 1000 for realistic young-galaxy spec¬ 
tra. Combined with the J 21 probability distribution by Di¬ 
jkstra et al. (2014), this gives usmbh ~ 10“^'’ Mpc“®, in 
rough agreement with the observed values. However, the pre¬ 
dicted nsMBH is still uncertain due to our poor knowledge of 
JJj from insufficient modelling of the physical/chemical pro¬ 
cesses in primordial gas clouds (see, e.g., Glover 2015a,b), 
and of the nature of the high-redshift sources responsible for 
this radiation. 

In the low-density primordial gas (< 10® cm“®), H 2 , 
the main coolant in the low temperatures regime, is formed 
either via the H~ channel or via the channel. The former 
begins with the formation of H“ by radiative attachment, 

H-be^H--b7, (1) 

followed by H 2 formation by associative detachment 

H^-bH-^Hz-be. (2) 

Similarly, the Hj channel starts with Hj formation by ra¬ 
diative association 

H + H+^H++7, (3) 

followed by the charge transfer reaction 

Hj + H Hz -b H+ . (4) 

In each channel, the second step proceeds much faster than 
the first one and its rate determines the amount of Hz for¬ 
mation. Since the rate coefficient for reaction (1) ~ 

6 X 10“^®cm®s“^ at 8000 K; Galli & Palla 1998, hereafter 
GP98) is about one order of magnitude larger than that for 

reaction (3) ~ 1-5 x 10“^®cm®s“^ at 8000 K; GP98), 

the H“ channel is usually much more efficient than that in¬ 
volving H^. However, under a strong radiation field, both 
channels become irrelevant due to the destruction of the in¬ 
termediaries by the H~ photodetachment 

H”-b 7 ^H + e, (5) 

and H^ photodissociation 

H++7^H + H+. (6) 

Since the binding energy of H“ (0.76 eV) is smaller than 
that of H^ (2.7eV for the ground state), H“ is photode- 
stroyed more easily than H^ and reaction (4) can become 
the dominant source of Hz in the primordial gas. 

Here, it should be noted that the H^ photodissocia¬ 
tion rate depends sensitively on the internal level popula¬ 
tion, because Hj in high vibrational levels is much more 
vulnerable to photodissociation. As an example, the bind¬ 
ing energy of H^ in vibrational levels v = 6 and 18 is 
only 1.2 and 0.003 eV, respectively. Thus, the level popu¬ 
lation of HJ is determined by a complex interplay of re¬ 
actions, involving formation (reaction 3), collisional disso¬ 
ciation (reaction 4), photodissociation (reaction 6), colli¬ 
sional excitation/de-excitation and radiative de-excitation. 
Therefore in order to accurately implement the effects of the 

J 21 = Jv{hu = 12.4eV)/10“^^ ergs“^ cm“^ Hz“^ sr“^ to quan¬ 
tify the intensity of external radiation. 


HJ channel, a detailed calculation of its level population is 
needed. 

In the post-recombination era, Hz formation is primar¬ 
ily controlled by the Hj channel for redshift 2 >100, owing 
to the suppression of H“ by the cosmic microwave back¬ 
ground (GMB) radiation (see, e.g., GP98). However, the 
efficiency of the Hj channel depends strongly on its level 
population: if HJ is in the ground state, the H^ channel is 
so efficient as to produce most of the Hz molecules at a level 
/iTz ~ 10“^. On the other hand, under local thermodynamic 
equilibrium (LTE), the H^ channel makes only a minor con¬ 
tribution to Hz, whose final abundance remains limited to 
/iTz ~ 10“®. Recently, Hirata & Padmanabhan (2006), Gop- 
pola et al. (2011) (hereafter Cll) and Longo et al. (2011) 
have studied the chemistry of the early Universe comput¬ 
ing the population of the vibrational levels of Hz and Hj 
following a state-to-state reaction kinetics, and found that 
HJ forms by reaction (3) preferentially in excited states. As 
a result, Hj channel does not dominate the formation of 
Hz, and the final Hz abundance is limited by H~ at a level 
fH, ~ 10 - 6 . 

Although the efficiency of the Hj channel depends 
strongly on the H^ level population, the LTE rate (Stancil 
1994; Mihajlov et al. 2007) has been widely used in studying 
primordial gas clouds without a real justification. An excep¬ 
tion is Glover (2015a) who calculated J 21 for black-body- 
type spectra with temperatures T^d ~ 10^ K and 10® K un¬ 
der two assumptions for the Hj level population: (i) all the 
H 2 + is in the vibrational ground state, and (ii) all levels are 
in LTE. Although the difference in Jzj in the two cases is 
not significant, it is not yet clear in what circumstances and 
to what extent the H^ channel can affect the evolution of 
primordial gas clouds. 

In this paper, we compute the thermal and chemical 
evolution of primordial gas clouds under a strong external 
radiation field with J 21 around Jzj by computing the H^ 
vibrational level population, and assess the effect of Hz for¬ 
mation via the H^ channel. We also determine the critical 
intensity for supermassive star formation and examine its 
dependence on the HJ level population by comparing our 
non-LTE results to those obtained in the LTE or ground 
state approximations. 

The paper is organized as follows. In Sec. 2, we describe 
our model for collapsing primordial gas clouds. In Sec. 3, we 
present the result of our calculation. The implications and 
conclusions are described Sec. 4. 

2 MODEL 

2.1 Basic Equations 

To follow the gravitational collapse of primordial gas clouds, 
we use the one-zone model described in SOI14 (see also 
Omukai 2001)), updated as follows: the vibrational level 
population of HJ is resolved following Gil; a part of the 
chemical network is updated following Glover (2015a,b). In 
our model, we compute physical quantities in the homoge¬ 
neous central part of the self-similar solution of collapsing 
clouds (Pension 1969; Larson 1969; Yahil 1983). The qual¬ 
itative validity of the one-zone model has been confirmed 
recently by three-dimensional hydrodynamical simulations 
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(Shang et al. 2010; Latif et al. 2014). Thanks to the sub¬ 
stantial simplification of gas dynamics, we can focus on the 
thermo-chemical processes in detail. 

According to the one-zone model, the evolution of the 
gas density p is modeled as 


^ = JL 

dt tg 


( 7 ) 


where tg = ^3tv/32G{p + pdm) is the free-fall time, G the 
gravitational constant and pdm the dark matter density, 
which is assumed to evolve following the solution of the 
spherical top-hat collapse model until it reaches the virial 
density. The radius of the core is approximately given by 
the Jeans length Aj = y/RjG^g^R/Gpp/trip, where rrip is the 
proton mass, p the mean molecular weight and k the Boltz¬ 
mann constant. The evolution of the gas temperature Tgas 
is determined by the energy equation. 


dCkin p . dCint 

~dr ~ “ ~dr 


( 8 ) 


where eidn = 3fcTga,s/2pmp is the kinetic energy of the gas 
per unit mass, Tc = (pfcTgas/Atmp) (d/dt) (1/p) the com- 
pressional heating rate, Arad the net cooling/heating rate 
due to radiative processes and eint the internal energy of 
the gas per unit mass, including both chemical energy and 
internal molecular energy. In the density and temperature 
range of our interest, the dominant cooling processes are H 2 
(Glover & Abel 2008), updated according to Glover (2015a), 
and Ly-a emission (Anninos et al. 1997). Other radiative 
reactions or the time variation of eint hardly affect the evo¬ 
lution of Tgas- 

We model the external radiation field as a diluted black- 
body with cut-off energy at 13.6 eV, specifying its tempera¬ 
ture Trad and intensity J 21 at 12.4 eV. We limit our analysis 
to black-body spectra with various T^ad although realistic 
spectra of galaxies are more complex and not well repre¬ 
sented by a black-body spectrum (SOI14, Agarwal & Khoch- 
far 2015). SOI14 claimed that the ratio of the photode¬ 
tachment rate to the H 2 photodissociation rate, fcpd /k^g, 
is a key parameter characterizing the hardness of the spec¬ 
tra, and that typical young and/or metal-poor galaxies have 
hard spectra with kfd /^pd 10^1 corresponding to black- 
body spectra with T^ad ^ 2 x 10"^ K. Nonetheless, we explore 
a broad range of black-body temperatures, Trad > 5 x 10® K, 
to mimic the evolution in a variety of environments possibly 
present in the early Universe. 


2.2 Chemistry 

In this section, we review our vibrationally-resolved chemical 
network based on SOI14 and Cll. We consider 29 chemical 
species (H, H+, H“, H 2 , H^(u) with u = 0, 1 ..., 18, He, 
He”*", He^^, HeH^, and e), and we compute the evolu¬ 
tion of the fractional abundance of species i, y{i) = n{i)/n^, 
where n(i) is the number density of species i and nn the 
number density of hydrogen nuclei. We do not consider deu¬ 
terium chemistry because it should only affect the evolution 
of the gas with Tgas below a few hundred kelvin by HD 
cooling, a temperature regime that is never attained in our 
calculations (see, e.g., Nagakura & Omukai 2005; McGreer 
& Bryan 2008; Wolcott-Green & Haiman 2011; Nakauchi 
et al. 2014). 



1000 10000 
Trad [K] 


Figure 1. Photodissociation rate of from the ground vibra¬ 
tional state (dashed red curve) and assuming LTE level popula¬ 
tions with = 1000, 3000 and 8000 K (solid red curves, bottom 
to top). The dashed green curve shows the H“ photodetachment 
rate. The radiation spectrum is a diluted black-body normalized 
as J 21 = 1. 


The processes involving HJ are the following. The vi- 
brationally resolved reactions in our code are H^ formation 
by radiative association (reaction 3), Hj photodissociation 
(reaction 6), dissociation by charge transfer (reaction 4), 
vibrational excitation/de-excitation by collision with H and 
vibrational de-excitation by spontaneous emission. We do 
not resolve rotational levels of H^, because the dependence 
of the efficiency of photodissociation on rotational levels 
is weaker than on vibrational levels (Dunn 1968; Babb 2014). 
In addition, complete rovibrational state-resolved data are 
not available in the literature. Note that HJ formation via 
HeH"*' hardly affects the evolution of primordial clouds (see, 
however, GP98 or Cll for its consequences on the chemistry 
of the early Universe). 

The state-resolved data adopted in our chemical net¬ 
work are the following. Babb (2014) provides rovibrationally 
resolved data for radiative association and HJ photodisso¬ 
ciation, which are summed up with respect to rotational 
levels to obtain the vibrationally resolved rate coefficients. 
To obtain vibrationally resolved HJ photodissociation rate 
we assume rotational level population as explained in Ap- 

pendix A. Fig. 1 shows the Hj photodissociation rates k_^^ 
for a black-body radiation field and J 21 = 1, obtained as¬ 
suming (rovibrational) LTE population with Trad = 1000, 
3000 and 8000 K, and all H^ in the ground state. Fig. 1 
clearly shows that Hj becomes more easily photodissoci- 
ated as the HJ vibrational state becomes more excited. 
Krstic’s database (Krstic et al. 2002; Krstic 2005; Krstic 
& Schultz 1999; Krstic et al. 2002; Krstic & Janev 2003) 
provides vibrationally resolved data for charge transfer and 
vibrational excitation/de-excitation by collision with H. The 
rovibrationally resolved de-excitation rates of H^ by spon¬ 
taneous emission are taken from Posen et al. (1983). They 
are summed up with respect to rotational levels to obtain 
the vibrationally resolved rates. To see the dependence of 
the evolution of the clouds on the assumed H^ level popula¬ 
tion, in addition to the standard state-resolved runs (here¬ 
after non-LTE model ), we also perform runs using the (rovi¬ 
brational) LTE H^ photodissociation rate (hereafter LTE 
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model) and runs using the photodissociation rate for Hj 
in the (rovibrational) ground state (hereafter ground-state 
model). 

We review here our treatment of H“ photodetachment 
and indirect H 2 photodissociation (the so-called Solomon 
process). Fig. 1 shows the H“ photodetachment rate coef¬ 
ficient for a black-body radiation field normalized to 
J 21 = 1 as a function of Trad. Here the cross section by John 
(1988) is used. As seen in Fig. 1, the H“ photodetachment is 
generally more effective than H J photodissociation although 
the latter depends on the level population. For the H 2 
photodissociation rate coefficient, we adopt the formula by 
Wolcott-Green et al. (2011), 

= 1.4 X 10-^Vsh J 21 s-^ (9) 

where /ah is the self-shielding factor (for the explicit form, 
see eq. 10 of Wolcott-Green et al. 2011). For the radi¬ 
ation spectra considered in Fig. 5, Eq. (9) gives fcpj = 
1.4 X 10“^^ s“^ without self-shielding (/ah = 1). The fact 
that fcpj and are decreasing functions of Trad can be 
understood as the photodissociation of Hj and the pho¬ 
todetachment of H“ are relatively more effective for softer 
spectra, compared to the photodissociation of H 2 . Note that 
the exact form of /ah is not settled in the literature (see, e.g., 
Draine & Bertoldi 1996; Wolcott-Green et al. 2011; Richings 
et al. 2014; Hartwig et al. 2015), which may introduce an un¬ 
certainty in the value of by a factor of a few, as pointed 
out in SOI14 and Glover (2015b). Similarly, evaluating fcpj 
with the intensity at a single frequency hi/ = 12.4 eV, may 
introduce an uncertainty in the value of J 21 by a factor of 
a few if the spectrum changes substantially within the LW 
bands, as mentioned in SOI14. However, since the above un¬ 
certainties should not affect the conclusion of this paper, we 
leave these issues for future studies. 

We updated the chemical network of SOI14, following 
Glover (2015a,b). We include dissociative tunneling of H 2 
to an unbound state (Martin et al. 1996) and H collisional 
ionization with H (Lenzuni et al. 1991; Gealy & van Zyl 
1987) and He (Lenzuni et al. 1991; van Zyl et al. 1981). We 
replace the direct H 2 collisional dissociation rate with the 
fit given in Martin et al. (1996). 


2.3 Initial Conditions 

Following Omukai et al. (2008) and SOI14, we take as ini¬ 
tial conditions of the calculation the following values cor¬ 
responding to the turnaround redshift 2 = 16: Tgas = 
21 K, riH = 4.5 X 10“®cm“®, y{e) = 3.7 x 10“'^, i/(He) = 

8.3 X 10~^, J/(H 2 ) = 2 X 10“® and y{i) = 0 for the other 
species. Note that the initial chemical composition hardly 
affects the thermal evolution of the clouds: at the beginning 
of the evolution well within the initial adiabatic contraction 
phase, most of pre-existing H 2 is photodissociated under a 
strong LW radiation field with J 21 around JJ); y(e) settles to 
the value determined by the condition that the recombina¬ 
tion time scale is about the same as the dynamical (free-fall) 
one; and the abundances of the other species reach the val¬ 
ues in chemical equilibrium (Omukai 2001). 



Ujj [cm-3] 

Figure 2. Temperature evolution in a collapsing primordial gas 
cloud irradiated by a diluted black-body radiation with Tad = 
6000 K and J 21 = 0.05 as a function of the cloud density: non- 
LTE model (solid red), ground-state model (dashed blue) and 
LTE (dotted green). The vertical lines (long-dashed black) demar¬ 
cate boundaries between the low, intermediate and high-density 
regimes. 

3 RESULTS 

In the following, we show how the temperature evolution of 
the cloud is modified by the assumed H^ vibrational level 
population. Then, we calculate the critical intensity of the 
radiation field J 21 as a function of the black-body tempera¬ 
ture and assess the importance of H 2 formation via the Hj 
channel on the evolution of the primordial clouds. 

3.1 Thermal evolution of the cloud 

The run of the temperature, chemical abundance and H^ 
level populations of a cloud irradiated by a black-body ra¬ 
diation field with Tad = 6000 K and J 21 = 0.05 as a func¬ 
tion of the H number density are shown Figs. 2, 3 and 4, 
respectively. In Figs. 2 and 3 we also plot for comparison 
purposes the results for the ground-state and LTE models, 
along with the fiducial results for the non-LTE model. In 
Fig. 4, instead, we also plot the LTE populations with the 
temperature given by Tgas in the non-LTE model (i.e., the 
solid red curve in Fig.2) along with the results for the non- 
LTE model, because the difference in Tgas in the non-LTE 
and LTE models makes it difficult to extract the effect of 
non-LTE chemistry by directly comparing the populations 
in these models. The results for the LTE model agree well 
with those in SOI14, the differences being due to the up¬ 
dated chemical network. Below, we describe the time (or, 
equivalently, density) evolution of the cloud with particular 
attention to the Hj vibrational level population. 

3.1.1 Low-density Regime: na < 1 cm“® 

In this density range, the temperature is still very low, 
Tgas < 10® K, although it rises with density by adiabatic 
compression, as seen in Fig. 2. Thus, almost all H^ is in the 
ground vibrational state if LTE is assumed (Fig. 4). How¬ 
ever, in the non-LTE model a finite amount of H^ is in 
higher vibrational levels because the HJ molecules tend to 
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Figure 3. Abundances of H 2 (red), (blue) and H“ (xlO“^; 
green) for the same cases as Fig. 2: non-LTE model (solid), 
ground-state model (dashed), and LTE (dotted) models. 

be formed by radiative association in excited states (reac¬ 
tion 4; Ramaker & Peek (1976); Babb (2014)). This effect is 
more pronounced at low temperatures where the v = 15 level 
is more populated than the v = 5 and 10 levels. Note that 
not all the settles to the ground state even though the 
collisional excitation/de-excitation rate is much smaller than 
the radiative de-excitation rate in this low-density regime. 
This is because not only radiative de-excitation, but also 
radiative association and HJ photodissociation contribute 
to the level population. In the non-LTE model, since 
a fraction of Hj is in high vibrational levels, the larger 
photodissociation reaction rate results in a smaller abun¬ 
dance than in the LTE or ground-state models (Fig. 3). As 
a result, the amount of H 2 formed is smaller in the non-LTE 
case. However, in this density range, the H 2 abundance is 
very low, //fj ~ 10“®, in all models and does not affect the 
thermal evolution. 

3.1.2 Intermediate-density Regime: nn ~ 1 — 10^ cm“® 

In this density regime, the evolutionary paths bifurcate into 
those of atomic and H 2 cooling depending on the assumed 
level population (Fig. 2). As seen in Fig. 4, at high tem¬ 
perature (-^ 4000 K), in the LTE approximation the excited 
levels are somewhat populated around ~ 10cm“^. However, 
in the non-LTE treatment only a smaller amount of H^ re- 
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Figure 4. Populations of in the n = 0, 5,10, and 15 vibra¬ 
tional levels (black, red, green and blue, respectively) in the non- 
LTE model (solid) for the same case as Fig. 2. The LTE popu¬ 
lations (dot dashed) with the temperature given by Tgas in the 
non-LTE model, i.e., the solid red curve in Fig.2, are also plotted. 

sides in excited states, because of the higher photodissocia¬ 
tion rate from those levels. In LTE the fractional abundance 
of excited levels, and thus the HJ photodissociation rate are 
overestimated. This results in smaller amounts of both Hj 
and H 2 , compared to the non-LTE model (Fig. 3). Due to 
the lower H 2 fraction of the LTE case, the cloud collapses 
along the atomic cooling track, whereas in non-LTE it fol¬ 
lows the molecular hydrogen cooling path (Fig. 2). Finally, 
in the ground-state model, the smaller H^ photodissociation 
rate and the consequent larger amount of H 2 (Fig. 3) cause 
an earlier onset of H 2 cooling relative to the non-LTE case 

(Fig. 2). 

3.1.3 High-density regime: nn > 10^ cm~® 

The different evolutionary paths found in the previous 
regime are not modihed at higher densities by the details of 
the various chemical processes. Considering the large tem¬ 
perature difference between LTE and the other models, it is 
not possible to distinguish between the effects of the level 
population and those related to the chemical abundances 
shown in Fig. 3. In this regime, unlike the LTE population, 
the excited levels are more populated than the ground state, 
as seen in Fig. 4. The reason is as follows: due to the high 
density, the H^ photodissociation is negligible compared to 
collisional processes, such as radiative association, charge 
transfer and collisional excitation/de-excitation. Since not 
only collisional excitation/de-excitation, but also radiative 
association and charge transfer continue to play important 
roles, the level population does not converge to the LTE 
value. This effect can be considered as a pumping mechanism 
that depletes the ground level and allows for the occupation 
of higher excited states. 

3.2 The critical LW intensity J 21 

Fig. 5 shows the critical LW intensity J 21 for the black-body 
spectra with temperature Trad. We find J 21 by using the 
bisection method, in which we examine whether the given 
J 21 is high enough for the cloud to collapse along the atomic- 
cooling track by totally suppressing H 2 cooling (for more 
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Figure 5. Critical intensity for the direct collapse J 21 ^ func¬ 

tion of the black-body temperature Trad of fh® irradiation radia¬ 
tion. In addition to the value in the non-LTE model (solid red), 
those in the ground-state (dashed blue) and LTE (dotted green) 
models are plotted. 
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Figure 6. The threshold value (dashed green), given by 
Eq. (14), above which the channel is more effective than 
the H~ channel at njj = 10^ cm~^. The critical intensity JjJ for 
the direct collapse shown in Fig. 5 (solid red) is overplotted. The 
explanation for the domains A-D is given in the text. 


details, see, e.g., SOI14 and references therein). In addition 
to our fiducial result for the non-LTE model, those for the 
ground-state and LTE models are also plotted in Fig. 5 for 
comparison. 

Two features in Fig. 5 should be noted. First, the slope 
of J 21 changes discontinuously around Trad ~ 7000 K, ac¬ 
companying the shift of dominant H 2 -formation reaction de¬ 
termining J 21 from the Hj channel (Trad 7000 K) to the 
H“ channel (Trad ^ 7000 K). We will discuss the condition 
for the HJ channel to be dominant in Sec. 3.3. Second, in 
the range Trad ^ 7000 K, J 21 depends on the assumed 
vibrational populations. Therefore, it is necessary to take 
into account the non-LTE level population to obtain 
J 21 correctly because the efficiency of the Hj channel, the 
dominant H 2 formation channel in this range, is sensitive 
to its vibrational level population. Compared to the correct 
J 21 by the non-LTE modelling, those for the ground-state 
and LTE models are smaller or larger, respectively. This is 
because in the ground-state model (the LTE model), com¬ 
pared to the non-LTE model, the photodissociation rate 
is lower (higher) around the density where the bifurcation 
into the two evolutionary paths occurs, and thus more (less) 
J 21 is needed to suppress the H 2 cooling. 

On the other hand, J 21 for Trad ^ 7000 K is insensi¬ 
tive to the HJ level population because the H“ channel 
rather than the Hj channel is important in this range to 
determine J 21 . Recall that typical high- 2 : galaxies have hard 
spectra corresponding to a black-body with Trad 2 x 10"^ K 
(SOI14), as mentioned in Sec. 2.1. Thus, our adoption of 
the LTE photodissociation rate in SOI14 to obtain 
JS for realistic spectra is justified. We plot J 21 only for 
Trad < 10"^ K in Fig. 5, because J 21 does not depend on 
the vibrational populations for Trad > 10“^ K. For J 21 in the 
range Tad > 10^ K, we refer readers to, e.g., Fig.3 of SOI14. 


3.3 Condition for the H^-channel dominance 

Here, we discuss the condition that the H^-channel domi¬ 
nates over the H~ channel. We introduce Jji, the threshold 
value of J 21 above which the channel is more effective 
in H 2 formation than the channel. 


To begin with, we estimate the H 2 formation rate via 
the H~ channel. The H“ fraction is determined by the bal¬ 
ance between the radiative attachment (reaction 1), asso¬ 
ciative detachment (reaction 2) and photodetachment (re¬ 
action 5) as: 


?/(H-) = 


kfs. y(e) 


fcST+C'MH) 

The H 2 formation rate via the H~ channel is given by 

dySr 


( 10 ) 


dt 


= Cd 2 /(H-)n(H) 

^ fcra y(e)n(H) 

“ l + fc?dVKH)fcSi" 


( 11 ) 


where we have used Eq. (10) in the second fine. Similarly, 
we estimate the H 2 formation rate via the Hj channel. The 
fraction is determined by the balance between radia¬ 
tive association (reaction 3), charge transfer (reaction 4) and 
photodissociation (reaction 6) as: 


y(H2+) = 


1/(H+) 

C/+CV«(H) 


( 12 ) 


The H 2 formation rate via the Hj channel, i.e., the charge 
transfer rate, is given by 

J H+ 

dfe, . H+ 


dt 


= fc,? y(H+) n(H) 

_ fcfa^ y(H+) n(H) 

l + k^I/(n(H)k^J 


(13) 


where we have used Eq. (12) in the second line. By equating 
Eqs. (11) and (13), we obtain 


■^21 — 




- 1 


k^A nn 


kpd ^21=1 


(14) 

where we have used fcpj = fcpj |j 2 i=i >^ 21 , ’T'(H) ~ uh 
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JJ + 

and y(e) ~ and neglected for simplicity. The 

effect of photodissociation would substantially change 
the expression given by Eq. (14) if the LTE rate was used. 
However, the actual rate of photodissociation just before 
the atomic- and H 2 -cooling tracks bifurcate is lower than the 
LTE one, as shown in Sec. 3.1.2. Thus, the correction due 
to HJ photodissociation makes Jji given by Eq. (14) larger 
at most by a factor of a few and can be neglected in the 
order-of-magnitude estimate here. 

To evaluate Jf}, we consider the physical condition 
of the gas just before the atomic- and H 2 -cooling tracks 
bifurcate (see, e.g., Omukai (2001) or SOI14) and take 
Tgas = 8000 K and uh = 10^ cm“®. The density nu at bifur¬ 
cation depends on the spectra and it is closer to 10^ cm“® 
for soft spectra with Trad ^ 10^ K, which is of our in¬ 
terest in this paper, although it is about 10^ cm~® for 
Trad 10® K. The reaction rate coefficients at 8000 K are 
given by « 6 x 10~^® cm® s“®, « 9 x 10~^® cm® s“® 

and k^^ ~ 1.5 x 10“®® cm® s“® (GP98) and the H“ pho¬ 
todetachment rate is given by Fig. 1. By substituting the 
above reaction rate coefficients and photodetachment rate 
into Eq. (14), we obtain J 21 , as shown in Fig. 6. Also shown 
is J 2 I obtained in Sec. 3.2 for comparison. Depending on the 
value of J 21 relative to JJi and the way in which a cloud 
evolves falls within one of the following four cases (Fig. 6): 

A: J 21 > J 21 and J 21 > Jf}. The cloud collapses along the 
atomic-cooling track, since H 2 , formed mainly via the Hj 
channel, is not enough to cool the gas. 

B: J 21 > J 21 and J 21 < .Rf The cloud collapses along the 
atomic-cooling track, since H 2 , formed mainly via the H“ 
channel, is not enough to cool the gas. 

C: J 21 < J 21 and J 21 > The cloud collapses along the 
H 2 -cooling track due to H 2 formation mainly via the 
channel. 

D: J 21 < J 21 and J 21 < .Rf The cloud collapses along the 
H 2 -cooling track due to H 2 formation mainly via the H“ 
channel. 

In the soft-spectrum regime of T^ad ^ 7000 K, where Jji < 
J 21 , the cases A, C, and D can be realized depending on 
J 21 . In the case C, in particular, while the H“ channel is 
blocked by radiation, sufficient H 2 for cooling can still form 
via the channel. In this case, the proper account of the 
level population is indispensable since the H 2 formation 
rate via the Hj channel is sensitive to its level population, as 
shown by the example in Sec. 3.1. On the other hand, in the 
hard-spectrum regime of Trad ^ 7000 K, where J 21 < J|i, 
the cases A, B and D are allowed: there is no range of J 21 
satisfying Jj)® < J 21 < J 2 i- Namely, radiation fields strong 
enough to block the H“ channel always exceed the critical 
intensity to totally suppress H 2 formation by H 2 photodis¬ 
sociation. This means that the H 2 -cooling track cannot be 
realized by the channel for fyad ^ 7000 K. 


4 DISCUSSION AND CONCLUSIONS 

We have computed the thermal and chemical evolution of 
primordial gas clouds resolving the vibrational levels of H^, 
which enables us to properly implement H 2 formation via the 
HJ channel. The efficiency of HJ photodissociation, which 


suppresses the channel by destroying the intermediate 
product H^, is sensitive to the level population of the molec¬ 
ular ion. We have found that H 2 formation via the chan¬ 
nel becomes more effective in the non-LTE model than in 
LTE because more Hj is in the ground state and thus its 
photodissociation rate is smaller. 

As to the effects of the background radiation, we have 
found that H 2 formation via the channel becomes im¬ 
portant in clouds irradiated by strong radiation fields with 
soft spectra characterised by fyad fy 7000 K. In this case, the 
cloud thermal evolution is strongly affected by the level 
population, indicating the importance of non-LTE treat¬ 
ment of the HJ vibrational levels. On the other hand, under 
radiation fields with harder spectra and fyad ^ 7000 K, the 
HJ channel always falls short of the H“ channel in produc¬ 
ing H 2 molecules. 

We have derived the critical radiation intensity J 21 for 
the formation of supermassive stars by direct collapse and 
examined its dependence on the assumed level popu¬ 
lation. For soft spectra with fyad fy 7000 K, the critical 
intensity J^j is under(over)-estimated if the level popula¬ 
tion is assumed to be in LTE (if all the Hj is assumed to 
be in the ground state). On the other hand, the value of 
J 21 is independent of the assumed level population for 
harder spectra with fyad ^ 7000 K. Therefore, since typical 
high-redshift radiation sources, i.e., young and/or metal- 
poor galaxies, have rather hard spectra corresponding to 
black-bodies with fyad // 2 x 10^ K (SOI14), the LTE ap¬ 
proximation adopted in SOI14 in determining J 21 for real¬ 
istic spectra is justified. 

We note that primordial gas clouds may be exposed 
to much softer radiation fields. For example, galaxies with 
extremely strong Ly-a emission can be regarded as charac¬ 
terised by very soft spectra. Since the energy of Ly-a pho¬ 
tons {hv = 10.2 eV) is below the lower limit of the Lyman- 
Werner bands [hu = 11.2 — 13.6eV), such emission would 
only destroy H“ leaving H 2 unaffected. Also, since the ini¬ 
tial mass function in high-redshift galaxies is poorly con¬ 
strained, the infrared/optical cosmic background radiation 
originating from these sources can be soft with temperature 
Trad ~ 6000 K, corresponding to the typical stellar mass 
~ 0.7 Mq, as considered in Wolcott-Green & Haiman (2012). 

In non-standard cosmology primordial density fluctu¬ 
ations are expected to be enhanced at small scales, so 
that the first stellar objects can form at very early times 
with 2 : > 100, as considered in Hirano et al. (2015). Un¬ 
der the strong and very soft CMB radiation at that time 
(Trad ^ 300 K), H 2 formation via the channel can be 
the dominant path with the H“ channel blocked by radia¬ 
tion processes (including non-thermal photons that play an 
additional role in enhancing H“ photodetachment; Coppola 
et al. 2013). 

As a last example, let us consider the case where direct 
collapse via the atomic-cooling track is realized. In such a 
case, the growing protostar results in a very extended struc¬ 
ture with radius R ~ 30 AU (“super-giant” protostar) and 
low effective temperature fyad — 5000 K (Hosokawa et al. 
2012, 2013). In studying the physical conditions in the re¬ 
gion around the super-giant protostar, the non-LTE level 
population must be taken into account. 

In this paper, the collapse of primordial gas clouds has 
been calculated by resolving the HJ vibrational levels, to 
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our knowledge, for the first time. However, neither rota¬ 
tional levels nor H 2 rovibrational levels have been explicitly 
included. Although the rotationally-resolved kinetics 
is probably not important, the H 2 rovibrationally-resolved 
kinetics may have considerable effects. The H 2 rovibrational 
level population affects critical processes in the evolution 
of the primordial gas, such as photodissociation, collisional 
dissociation, and cooling by molecular hydrogen. We plan 
to include the full level-resolved kinetics of H 2 and in 
future studies. 
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Figure Al. Same as Fig. 5, but for three different assumptions 
on the rotational level population: (*) model using the inter¬ 
polation formula by Glover (2015a), {ii) model where all the Kf 
in each vibrational level is in its rotational ground state, and (in) 
model assuming the rotational LTE. Note that the vibrational 
levels are resolved in all the cases. 
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APPENDIX A: Uf ROTATIONAL LEVEL 
POPULATION 

Here, we describe our treatment of the HJ rotational level 
populations. Since comprehensive rotational-state-resolved 
data is not available in the literature, we use the interpo¬ 
lation formula by Glover (2015a). Specifically, we calculate 
the H^ photodissociation rate from each vibrational level by 
the formula k = ^ltb (fcn->o/fcLTB)“, where fcn-»o and ^ltb 
are the rates for the ground rotational state and the LTE 
rotational level population, respectively, a = (U-UH/ucr)”^ 
and Ucr is the critical density for the LTE populations, which 
is ~ 10® cm“® for Tgaa = 8000K (Glover 2015a). 

To estimate the error in our modeling of the H^ ro¬ 
tational level population, we compare the critical intensity 
J 21 obtained with three different modeling of rotational level 
population: (i) the fiducial model based on Glover (2015a) 
(j-interp. model); (ii) the model assuming that all the Hj in 
each vibrational level resides in the ground rotational state 
(j'-ground state model); and (hi) the model assuming the 
LTE rotational level population in each vibrational level (j- 
LTE model). Fig. Al shows that the values of J 2 I in the j- 
interp. and j-ground state models are almost identical while 
in the j-LTE model it is slightly smaller. This difference is, 
however, much smaller than that among the vibrationally 
non-LTE and LTE models shown in Fig. 5, indicating that 
the cloud thermal evolution does not depend so much on ro¬ 
tational level population as the vibrational level population. 
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